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Abstract. In this paper explore the relation between covariant and canon- 
ical approaches to quantum gravity and BF theory. We will focus on the 
dynamical triangulation and spin- foam models, which have in common that 
they can be defined in terms of sums over space-time triangulations. Our aim 
is to show how we can recover these covariant models from a canonical frame- 
work by providing two regularisations of the projector onto the kernel of the 
Hamiltonian constraint. This link is important for the understanding of the 
dynamics of quantum gravity. In particular, we will see how in the simplest 
dynamical triangulations model we can recover the Hamiltonian constraint via 
our definition of the projector. Our discussion of spin-foam models will show 
how the elementary spin-network moves in loop quantum gravity, which were 
originally assumed to describe the Hamiltonian constraint action, are in fact 
related to the time-evolution generated by the constraint. We also show that 
the Immirzi parameter is important for the understanding of a continuum limit 
of the theory. 



General relativity can be quantised by canonical or covariant methods depending 
on whether one regards space or space-time as more fundamental. Because of the 
incomplete nature of the subject the relation between the two approaches is not 
entirely clear, even though one expects that in appropriate circumstances they 
should be equivalent. The broad aim of this paper is to contribute to this discussion. 

General relativity is a completely constrained system and hence quantum grav- 
ity distinguishes itself from ordinary field theory in that the Hamiltonian vanishes 
weakly. This implies for a canonical formalism that there is no propagator, since 
there is no external time variable with respect to which it could be defined. In- 
stead the central object is the generalised projector onto the solution space of the 
Hamiltonian constraint H(x). Formally, this projector can be written as: 



where [dT(x)] is a measure on the space of scalar fields T(x) and E is a spatial slice 
of dimension n of the space-time manifold M. = E x M. This operator then defines 
the physical Hilbert space of the theory, from which we have to extract all physical 
predictions by the construction of suitable observables. 

This projector can be expressed as a phase space path integral in an analogous 
way as is done for the propagator in ordinary field theory. Hence, we expect this 
projector to be the link between the above canonical description and covariant 
approaches to the quantisation of gravity. 
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Indeed, the starting point for the covariant approach is the expression 

f92 

Z(gx,g 2 ) ■■= / [dg]e-^ 

J 91 

which describes a sum over all space-time geometries g on the cylinder S x [0, 1] 
with boundary geometries g\ and gi. Two promising approaches towards making 
sense of this expression are the dynamical triangulation and the spin-foam models. 
These have in common that they they replace the above sum over geometries by a 
sum over space-time triangulations. 

The aim of this paper will be to show that the projector P can be suitably 
regularised so that it describes a sum over space-time triangulations. More precisely, 
we will be discussing canonical theories of quantum gravity that have a Hilbert space 
of kinematical states A given by triangulations or dual graphs of the spatial slice E. 
We will provide two different regularisations of the projector such that its matrix 
elements (Ai |P| Ao) are given by a sum over triangulations of the cylinder Ex [0, 1], 
with boundary triangulations Ao and Ai. The amplitudes of the triangulations will 
be given by the dynamical triangulation and spin-foam models respectively. 

The idea will be to decompose the projector into a product of "small" unitary 
local proper time evolutions. Acting on an initial triangulation A each elementary 
evolution corresponds to a gluing of an (n+l)-dimcnsional simplex to the boundary. 
In this way we build up a sum over space-time triangulations with the correct 
amplitude. 

This link between canonical and covariant theories is important for the under- 
standing of the dynamics of these quantum gravity models and for providing a 
space-time picture in the canonical framework. Ultimately, one would like to relate 
a given Hamiltonian to the amplitude of a corresponding covariant model. As we 
will see this can be achieved in the simplest case of 1+1 dynamical triangulations. 
The discussion of spin foam models of quantum gravity and the closely related 
topological BF theories, will clarify their relation to canonical loop quantum grav- 
ity. In particular, we will see that elementary moves of spin-networks, which were 
thought to describe the action of the Hamiltonian constraint, arc in fact related 
to time-evolutions. A better understanding of this relation should help in deciding 
which of the current proposals for spin foam models provides a correct description 
of quantum gravity and also shed light on the nature of the Hamiltonian constraint 
of loop quantum gravity. 

This work builds on a previous attempt to derive a spin foam model from 3+1 
loop quantum gravity by Reisenberger and Rovelli pc| , which was developed further 
by Rovelli in b3] . Our approach will differ in the regularisation of the projector. 
We will see in this paper how this leads to many attractive features: 

• By providing two different regularisations we can describe both the spin-foam 
and the dynamical triangulation models. 

• Our expansions are closely related to the derivation of the path integral in 
quantum mechanics and field theory. 

• They provide a clear role for the [cff'(cc)] integration over scalar (lapse) func- 
tions in the projector and avoid infinities due to this integration. 

• They provide a very natural interpretation of the addition of space-time sim- 
plices as local proper time evolutions. 

• They avoid over counting of space-time triangulations. 

• We obtain an understanding of the continuum limit of spin-foam models. 
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• Most importantly one can show that in BF theories and 2+1 loop quantum 
gravity the addition of tetrahedra, when viewed as an operation on physical 
states, leaves these invariant. This suggests that these actions should be 
viewed as evolutions and not actions of the Hamiltonian constraint operator. 

• In addition, in 1+1 dynamical triangulation we can show explicitly how we 
can derive the correct Hamiltonian operator from the covariant model via our 
expansion of the projector. 

We begin our discussion by describing the formal relation between the projector 
and path integrals for quantum gravity. This will be made concrete in the dynam- 
ical triangulations framework focusing on 1+1 dimensions, where results can be 
stated most explicitly and simply. We then continue to generalise our results to 
higher dimensions before moving on to consider spin-foam models and loop quan- 
tum gravity. 

2. Canonical vs. Covariant 

In this section we describe the generalised projector and show how it is formally 
related to a sum over geometries. We will find that canonical and covariant quan- 
tum gravity are equivalent up to the choice of range in the integration over lapse 
functions. 

2.1. Canonical quantum gravity. The canonical description of general relativity 
requires that we restrict our attention to space-time manifolds M that allow a split 
into space and time, i.e. manifolds diffeomorphic to M. = £ x R. Here and in the 
following £ will denote the n-dimensional spatial slice, which we will take to be 
compact. 

General relativity is a completely constrained system, i.e. the Hamiltonian for the 
theory vanishes weakly This can be viewed as a consequence of diffeomorphism 
invariance and the absence of an external time parameter. Hence quantisation 
follows the Dirac prescription which requires that we quantise the unconstrained 
classical system and impose the constraints as quantum operators on an initial, 
kinematical, Hilbert space Hkin- Loosely speaking, the constraints implement the 
diffeomorphism invariance of the theory. Typically, we distinguish between the 
"momentum" constraints Hi, which generate diffeomorphisms of the spatial slice 
£ and the "Hamiltonian" constraint H, which is responsible for diffeomorphisms 
normal to the spatial slice and is hence associated with dynamics. If connection 
instead of metric variables are used in the coordinatisation of the phase space there 
will be additional "Gauss" constraints, enforcing invariance under internal gauge 
transformations corresponding to the rotation of tetrads. 

Physical states lie in the kernel of the constraint operators. If this space is non- 
trivial and can be endowed with an inner product, we obtain the final physical 
Hilbert space TL p hys ■ Note the absence of a Hamiltonian operator on TL p h V s , result- 
ing in the the lack of a notion of evolution or propagator. States are interpreted 
as entire histories, which makes the extraction of physical information from this 
frozen picture highly problematic. 

These conceptual matters will not concern us. Instead we need to face the 
following technical issue related to our problem of constructing a path integral for 
general relativity. Whenever lies in the continuous part of the spectrum of one of 
the constraint operators, the corresponding eigenvectors will not be normalisable. 
Hence the construction of the physical Hilbert space is non-trivial and proceeds by 
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the so-called refined algebraic quantisation (RAQ) procedure (see for a review 
and references). The idea is to look for distributional solutions to the constraints, 
i.e. functionals on the kinematical state spacejjj The solutions are precisely those 
functionals that are constant on the orbits of the unitary action generated by the 
constraints. Equivalently we could be view the orbits themselves as the generalised 
solutions. The key step is to define an inner product on these physical states, which 
is done by introducing a generalised projector P. This is a map from Jikin onto the 
space of functionals, which maps a vector to its orbit under the constraint action. 
The image of the map P will be referred to as the solution space of the constraints. 

Let us denote physical states by a tilde, i.e. ip := Pip, <p :— P<j> with tp, <f> G Hkin, 
then the physical inner product is given by: 

{ipA)phys ■= <f>[lp], (1) 

where the right hand side denotes the action of the functional <p on the kinematical 
state ip. In many cases P can be constructed by the method of group averaging. 
In this case we have: 

(ip,4>) P hys ■■= J dg{ip,U{g)(t>)kin ■■= (ip,P<t>)kin (2) 

where U (g) denotes the group action generated by the constraint under consider- 
ation and dg is a suitable measure on the group. We have also made use of the 
kinematical inner product (•, -)kin on Tikin- 

Completion of the solution space with respect to the physical inner product gives 
us the physical Hilbert space H p ) ly s- The operator P defined above is called the 
generalised projection operator since it reduces to an ordinary projection onto the 
kernel of the constraints if the corresponding group action is compact. But the 
power of the approach lies in the fact that it can be applied in other cases also, 
in fact whenever the above integral converges, group averaging defines the unique 
generalised projector ||. 

When applied to the momentum and Gauss constraints the above procedure is 
implemented simply, by demanding that states are invariant under the action of 
spatial diffeomorphisms and gauge transformations respectively. The Hamiltonian 
constraint, however, causes difficulties since it generates diffeomorphisms which are 
normal to the spatial slice and hence does not have a simple geometrical action on 
kinematical states. Indeed, the Hamiltonian constraint is related to evolution in 
coordinate time hence we expect the projector corresponding to H to be related to 
the propagator in field theory. This is what we explore in the following. 

In analogy to the group averaging procedure we can write down a formal ex- 
pression for the generalised projection operator corresponding to the Hamiltonian 
constraint operator H(x): 

P = J{dT(x)]e-^ T \ (3) 

where H[T] := J^d n xH(x)T(x). Here [dT(x)] denotes a measure on the space of 
scalar functions T(x), which we leave unspecified. The integral should be under- 
stood as an integral over the generators of the algebra generated by H(x). Since 
the exponential exp(— iH[T]) is linear in T(x), the above can also be viewed as a 



To be precise we need to consider a dense subspace of the kinematical Hilbert space. This 
technical issue will not be important for our formal discussion. 
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expression of a product of delta functions imposing the constraint Ilxes ^[^( x )]- 
Alternatively, we can understand the above in the following way. The expression 
(ip, e lH ^ 4>) kin defines the "propagator" in multi-fingered proper time T{x) between 
two kinematical states ip and <f>, which describes geometries on the boundaries of 
the cylinder [0, 1] x S. Since physical results should be independent of the arbitrary 
proper time separation between initial and final splices, we should integrate over all 
possible T(x). This then defines the above projector. A more in depth discussion 
is given in pfj| . We will now show in more detail how the projector and the path 
integral for quantum gravity are related. 

2.2. Covariant quantum gravity. The projector can be expressed as a phase 
space path integral as is demonstrated by Klauder |Q for general constrained 
systems. This phase space integral can also be derived from a sum over geometries 
approach to quantum gravity, as has been investigated in detail by Teitelboim [ pTj[ 
[3l| , forging the link between covariant and canonical frameworks. 

Let us denote the points in the phase space of general relativity as [p(x),q(x)], 
x G E and let us furthermore assume we have eigenstates \q) £ Ti-kin of the corre- 
sponding configuration operators. We can then express the physical inner product: 

{qx,qo)phy S := [ [dT(x)]{ qi \e- i6 ^\q ) (4) 



as a path integral in complete analogy with the derivation of the path integral 
in elementary quantum mechanics. More concretely, we imagine that the operator 
exp(— iH [T]) is evolving between a state describing an initial geometry at coordinate 
time and a final geometry at coordinate time 1. We then write: 

N 

e-^m = linjn^' T1 . (5) 

e_> fc=i 

with e = 1/N. Then using the standard techniques we obtain formally: 

(?l,«)W= fldT(x)} I"" 1 [ D p}[Dq}e^o dtj d^( P ( x ) t ^) t ~H( P , q )T( x) ) (g) 

J Jq 

The expression on the right is an integral over all paths in phase space 

t^[q{x) u p{x) t ] fe[o,i] (7) 

with fixed initial and final configurations go & n d qi- The measure [Dq] := ]^[ t _ [cJQ'(a;)] 
is a product of suitably normalised measures on the space of functions q(x) 

with a similar expression for [Dp] . 

This is the expression that is obtained if we start from a sum over histories 
approach to quantum gravity modulo ambiguities concerning the integration range 
in the [dT] integral. There one is interested in defining the quantity: 

f9l 

/ [dg]e- lS[sl (8) 
J go 

where S denotes the action for general relativity. The integral is over all space-time 
geometries^ g that interpolate between initial and final geometries go,9i on the 
boundary of the cylinder [0, 1] x E. By a geometry we mean an equivalence classes 
of metrics under diffeomorphisms. 



2 Here and in the following bold face characters will correspond to space-time quantities 
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This sum over geometries can be reexpressed as a phase space path integral. In 
metric variables the phase space coordinates are given by the metrics gi k (x) on the 
slice £ and their conjugate momenta p lk (x) . 

By choosing a foliation of the space-time cylinder we can rewrite the gravita- 
tional action in Hamiltonian form. The action is then a function of the paths 
t — > [gik{x, t),p lk (x, t)], where the pair [gik{x, t),p lk (x, t)] describes the geometry of 
the spatial slice {t} x £ 
1 

ik • /\r/ +\tt{~ ^'ik\ Aji(™ 4-\ TJ 



S[N,N%p,g] = / dt d n x(p lk g lk -N(x,t)H(g zk ,p lk )-N\x,t)H l (g lk ,p lk )) 
Jo Jt, 

(9) 

where N and iVj are the lapse and shift functions which encode the information on 
how the foliation of space-time is chosen. 

The constraints generate gauge transformations that relate physically indistin- 
guishable variables. Hence, to avoid multiple counting of equivalent paths in the 
integral (||) we need to fix a gauge. One of the simplest possible gauge choices, 
as described by Teitelboim is the proper time gauge N = 0, N z = 0. These 
conditions fix the available freedom up to diffcomorphisms of the initial or final 
slices. Hence one can rewrite the integral (|J) as an integral over paths in phase 
space with initial and final geometries g(x,0) and g(x, 1) kept fixed: 

/ [dT{x)] ]J[dp lk {x,t)}[dg lk {x,t)}e^o dt k d "^P' k 9^-Tix)H^ k y'')) (1Q) 
■* t=0 

where the integral is over the functions T(x) := N(x,t), which give the multi- 
fingered proper time separation between initial and final spatial slices. This integral 
has the same form as the one derived from the canonical projector in equation (||). 
Note however, that the integration measures need to be chosen in such a way that 
the above integral becomes independent on our choice of gauge, otherwise one 
needs to include ghosts terms in the above action as described by Teitelboim p9[ . 
Alternatively, by choosing the appropriate measures Klauder has derived the above 
directly, without the prior imposition of the gauge conditions Jl2| . We will not be 
concerned with these issues as the measures will be implicitly determined by the 
choice of models considered later, which provide concrete realisations of the formal 
path-integral expressions given in this section. 

2.3. Relating canonical and covariant approaches. We are interested in the 
range of integration of the [dT(x)] integral which we have left unspecified and 
which is responsible for a subtlety in the relation between canonical and covariant 
approaches. First we note that in order for the lapse function to arise from a well 
defined foliation it is necessary that T(x) ^ for all x. Since T(x) is a continuous 
function this implies that T(x) is either entirely positive or entirely negative. 

If the integral ([[o]) is to project onto the kernel of the Hamiltonian constraint 
and be equivalent to the expression @ then we need integrate over both positive 
and negative lapse functions, as is clear from the analogy between the projector 
and a product of delta functions. 

However, if we want the phase space integral ([To]) to be equivalent to the sum over 
geometries integral (||), then we should only integrate over half of the allowed T(x), 
say only those for which T(x) > 0. To see this we recall that that the function T(x) 
encodes how the space-time manifold M is foliated as [0, 1] x S. In equation ([l0|) we 
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are integrating over the variables [T(x), gik,p tk )- The function T(x) determines how 
the 3-metric on each slice {t} x E and their conjugate momenta are glued together 
to give a space-time metric on M.. However, a change in sign of the proper time 
T(x) simply corresponds to reversing the foliation of M. and hence [T(x), gik,P lk ] 
and [—T(x),gik,p lk ] will describe equivalent space-time metrics. 

This symmetry is not enforced by the constraints as they only identify config- 
urations that can be related by infinitesimal deformations. Hence the sum over 
geometries expression is not exactly equivalent to the phase space integral derived 
from the canonical picture due to a larger symmetry group in the sum over geome- 
tries case. This is reflected in the range of integration in the integral^ (|l0|) . 

This feature can be seen clearly in the discussion of the relativistic particle, which 
can be viewed as general relativity in 0+1 dimensions. If we choose the half infinite 
range of integration then the phase space path integral will define the Feynman 
propagator (an in depth discussion of the different Greens functions that can be 
obtained in this approach is given by Halliwell and Ortiz in (To) ) . This is the same 
result that is obtained if we try to define the sum over geometries integral, which is 
just a function of the length of the particle trajectory. One approach to do this is 
to use a dynamical triangulations approach in the Euclidean regime and then Wick 
rotate the result, which gives us precisely the Feynman propagator. 

Let us summarise. We have seen that P defined in equation (||) is an important 
operator for canonical quantum gravity. It has two interpretations depending on 
the T(x) integration range: 

1. For the full integration range we obtain the projector onto the kernel of the 
Hamiltonian constraint; 

2. For the half-infinite integration range P is related to the sum over geometries. 

In either case P can be described by a path integral. 

In the following we will show how this can be made concrete by discussing first 
the dynamical triangulations and second the spin-foam model of quantum gravity. 
The former gives a definition of a sum over geometries whereas the latter is a 
path-integral definition of the projector. We will see how these descriptions can 
be recovered from two respective regularisations of P defined on the appropriate 
kincmatical Hubert spaces. 

3. Dynamical triangulations 

The dynamical triangulations approach to quantum gravity attempts to regu- 
larise and define the path integral given by equation (g) (c.f. for details and 
references). The idea is to work with manifestly diffcomorphism invariant quanti- 
ties and thus avoid difficulties concerning gauge fixing. This is achieved by replacing 
the smooth space-time metric manifold by a simplicial manifold constructed from 
equilateral (n+l)-dimensional simplices of geodesic edge length a. Here and in the 
following we will be working with Riemannian metrics. A Lorentzian framework 
for dynamical triangulations exists j^j but for our purposes the Riemannian case 
will suffice to exemplify all the features of our construction. 

We are approximating the sum over geometries by a sum over triangulations of 
the space-time manifold. The weight of each triangulation is given by the Euclidean 



3 Teitelboim gives an argument why one might want to restrict the range of integration in the 
canonical case also Jjo| . The idea is that one should implement causality by demanding that one 
integrates only over configurations where the final slice is to the future of the initial one. 
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Regge action, which can be expressed solely as a function of the squared edge 
lengths. But note that in contrast to Regge calculus we are keeping the edge 
lengths fixed — triangulations differ only by the connectivity of the simplices. The 
challenge is to find a non-trivial continuum limit of this statistical model by taking 
a — > 0, and appropriately rescaling the bare Newton and cosmological constants. 
This works well in 1+1 dimensions, where one can recover results from standard 
Liouville theory. In higher dimensions evidence for a interesting continuum limit is 
still lacking. 

Our goal is to show how the sum over triangulations can be recovered from a 
canonical theory via an expansion of the projector defined in equation (|3|). We 
assume again that space-time is of the form M. — E x M. The kinematical Hilbert 
Hkin space for our theory is constructed from the free vector space generated by all 
equilateral triangulations A of E. By identifying the edges as geodesies of length a 
we can interpret states as piecewise-linear geometries. 

We will distinguish between auxiliary and kinematical states. Auxiliary states 
are actual embedded triangulations of E, whereas kinematical states are equivalence 
classes of auxiliary states under diffcomorphisms, i.e. we can think of them as 
non-embedded or abstract triangulations which are simplicial manifolds with the 
topology of E. 

An inner product is specified on the space of triangulations by demanding that 
incquivalent triangulations are orthogonal i.e.: 

(A, A') := *a,a' (11) 

Completion of the space of kinematical triangulations with respect to this inner 
product gives us Tihin- 

We are interested in the projector P acting on this space. From the discussion 
in the last section we know that P is related to the sum over geometries, which in 
the dynamical triangulation approach is given by 

Z Q (A 1 ,A ):=^e- 5[Al (12) 

which denotes the sum over all space-time triangulations ▲ with boundaries Ao and 
Ai weighted by the exponential of the Regge action S[A]. In particular, we expect 
that by restricting the integration range in the projector we should have: 

J [dT(x)]+(Ai \e-"^\ Ao) =Z o (Ai,A ) (13) 

where [dT(x)]+ denotes an integration over positive proper-time functions onlyQ 
Note that because we are doing a Riemannian calculation we are using real ampli- 
tudes. 

We now show how this can be demonstrated by choosing a suitable regularisation 
of the projector 



4 Note that this restriction implies that the expression on the left of equation ( Jl 3| ) does not 
define the physical inner product. 
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3.1. Regularisation of the projector I. Let us consider the projector with half- 
infinite integration range^ 



P = 



[dT(x)]+e- f d " xT ^^ 



dv / [dT(x)}+5 



d n xT(x) = v 



- f d n xT(x)H(x) 



(14) 
(15) 



Note that the dimension of v is naturally that of space-time volume. Hence we are 
splitting the sum over all geometries into a sum over geometries with fixed volume 
and then summing over all possible volumes. 

To regularise the functional integral, we allow the function T{x) to only take 
values in discrete steps of e. In other words T{x) = t(x)e with t(x) S N. We also 
write the spatial integrals as Riemann sums. To do this we assume that we have 
given some metric on E with respect to which we can define a square lattice with 
lattice spacing I. We will denote the set of vertices of this lattice by V. This gives 
us a regulated projector Pi >e . 



dv / [dT} e S 



**«£*( 



xev 



= £ / w 



s=0 



£*(*) 



.xev 



t{x) 



nn< 

xev k=o 



-d"H(x) 



where [dT] e denotes a measure on the positive discretised lapse functions. 
Let us set 



-trH(x) 



(16) 



(17) 



(18) 



One can see that the above delta function imposes the condition that the product 
rixev 111=0 Ui >e {x) has exactly s factors. The integral over the discretised lapse 
function then ensures that for fixed s we sum over all possible t(x) satisfying the 
constraint. More precisely we have: 



Pu = 



(19) 



s=o sev s 



where the second sum is over all possible ordered s-tuples x :— (x\, ...,x s ) € V s . 
The ordering arises in the replacement of the spatial integrals by Riemann sums. 
For example, in (l+l)d the ordering implies that x\ < xi < ... < x s . 

The continuum projector is recovered in the limits e, I — > 0, whereby we need to 
check that the results are independent of the lattice chosen in the regularisation. 

In the following we will see that this expansion allows a very natural interpreta- 
tion. The operator Ui i£ (x) will correspond to a local evolution of a state in proper 
time e, which will be identified with the addition of a space-time simplex to the 
boundary triangulation. In this way the above sum will correspond to a sum over 
triangulations of space-time with the number of tetrahedra in a particular triangu- 
lation given by s. In the next sections we will discuss this in detail in the context 
of 1+1 dimensions. 



5 Although this is strictly speaking not a projector we will, by abuse of notation, continue to 
refer to it as such. 



10 



MATTHIAS ARNSDORF 




Figure 1 . The evolution operator acts on a kinematical state by 
a adding a triangle in two possible ways. This creates a new kine- 
matical state which is depicted by the dashed line. 



3.2. 1+1 dimensions. In 1+1 dimensions the action for general relativity is pro- 
portional to the volume of space-time: S[g] = A J M y '— det g, where A denotes the 
cosmological constant. In the following we will be approximating the sum over 
geometries by a sum over triangulations of the cylinder S 1 x [0,1]. 

Space-time simplices are simply triangles with area proportional to a 2 . Hence 
the transition amplitude between two boundary states Ao and Ai is: 

Z a (A 1 ,A )-^e^ 2 ^ 

A 

where the sum is over space-time triangulations ▲ with the appropriate boundaries, 
and N A is the number of triangles in A. 

Auxiliary states correspond to triangulations of S 1 and are thus an embedded 
sequence of links. They can be characterised by the coordinates Vi of the vertices 
between the links, i.e. they are of the form \v) := \v\...v n ), where n is the number 
of links. Since the only diffcomorphism invariant information contained in the 
auxiliary states is the number of links, kinematical states can be though of an 
abstract chain of links denoted by |n). These are the sates spanning Ttkin- 

The aim is to find a to find a definition of Ui^{x) acting on Tlkm such that 
as I — ► we can express (n|P; jC |r7i) as a sum over triangulations with the correct 
weights. In other words we would like to deduce the appropriate quantisation of 
the evolution Ui ie (x) from the path integral form of the theory. To do this we will 
identify the action of the operator Ui t€ (x) with the stacking of triangles. 

There are two ways to add triangles to a sequence of links \n). We can glue 
one or two faces of the triangle to the boundary. This creates a new set of links 
|n+ 1) or \n — 1) respectively. Roughly, we will define Ui >e (x) to be the sum of both 
these local actions as shown in figure |l]. Repeated action of L\ e will then create 
a triangulation. An important aspect of this identification (which is unique to the 
dynamical triangulations framework) is that the parameter a in our model can be 
interpreted as a proper time distance. Hence it is natural to identify a with e, as 
we will do in the following. 

These ideas can be expressed more precisely in the dual picture, which is also 
more likely to be useful in generalisations to higher dimensions and spin-foams. In 
general, the dual of a (n + l)-dimensional triangulation is an (n + 2)-valent graph. 
The graph has a node in the centre of each simplex and one edge crossing each 
face as shown in figure |^. In one dimension the dual of a triangulation is just a 
sequence of nodes. Hence states that are dual to an auxiliary state \v±, ...,v n ) will 
be denoted by \x) := \xi, .., x n ) where Xi :— (wj + Vi + i)/2 now is the co-ordinate of 
the i'th node which corresponds to the centre of the i'th link. 

The operator Ui >a {x), being local, acts on auxiliary states. Let us introduce 
creation and annihilation operators a){x) and a(x) acting on nodes of auxiliary 
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Figure 2 . A 2-dimensional triangulation and its dual trivalent graph 




FIGURE 3. The actions of a t (x 2 ) (left fig.) and a(x 2 ) (right fig.) 
on the auxiliary state \xx,X2,Xa). The newly created states are 
depicted by the dashed lines and X1.5 — (xi + x%)j2. 



dual states: 

a+ / \ j l [ ifx 7^ Xi / nr\\ 

a'{x)\xi,...x n ) := < . xj-i+xj s . f _ (20) 

a(x)\xi,...,x n ) := I ? , \ •!■!!_ ( 21 ) 

This can be identified with the gluing of triangles as follows. The creation 
operator (v{xi) corresponds to the gluing of one edge of a triangle to a boundary 
link centred at x\. This gives a new boundary triangulation, which is identified 
with an auxiliary state by assigning the the node corresponding to right free edge 
of the added triangle with the coordinate Xi and the left edge with (xi-i + Xi)/2 
as depicted on the left of figure ||. 

The annihilation operator a(xi) acts by gluing two edges of a triangle to the 
boundary links Xi and a?i_i. This gives a new state with one less link than the 
original one. The free edge of the added triangle is identified with the node Xi-i 
as can be seen on the right of figure |[ We use these operators to define Ui ia (x): 

U lta {x)\x) := i e - Q2 V(.x) + a(x))\x) (22) 

Multiple actions of U^ a {x) will thus give rise to a sum of auxiliary states, each 
of which can be identified with an abstract triangulation. The projector is given 
by a sum of ordered products Ui, a (xi) ■ ■ ■ Ui. a (x s ) as shown in equation ([l9|). In 
the (l+l)d case the ordering means that x± < X2 < •■■ < x s . This implies that an 
ordered product of evolution operators will construct triangles from right to left (if 
the coordinate axis points left to right). 

The action of Ui, a {x) on a state \x) will be whenever x is not the coordinate of 
a node in \x). It follows that whenever the action of a product Ui, a {xi) ■ ■ ■ Ui, a (xs) 
is non-zero it will act by constructing a sum of triangulations, each containing s 
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triangles. Hence Pi >a corresponds to a sum over triangulations, where the parameter 
s in equation ( |l9| ) determines the number of triangles in a triangulation. 

As I — > the separation between the allowed values for x £ S 1 tends to 0. Hence, 
in this limit, any triangulation that can be constructed by an ordered sequence of 
evolution operators will be present in the sum over triangulations described by Pi a 

We now complete the definition of (n|Pj, a |m) as follows. We choose an auxiliary 
state |m) which is a representative in the class |m), i.e. any auxiliary state that has 
m links. The action of P/ l0 on this state will produce a sum of auxiliary states: 

Pl,a\m) =J2 e ~ a2XNli Wi)> ( 23 ) 

i 

where each of these states is the (future) boundary of an abstract triangulation 
k; containing A Ai triangles as described above. This gives: 

(n\P lia \m) =^e- Q2AAr ^(n,f l ) (24) 

i 

The inner product (n, x%) is given by the dual action of \n) on the auxiliary state 
\xi). This is 1 if if the number of links in \Si) is n and otherwise. 

In this way (n\P^ a \m) is given by a sum over triangulations that have two bound- 
aries, with n and m links respectively. The amplitude for each triangulation is just 
e -a \N k ^ w } iere jy A j s the- number of triangles in the triangulation as required. 

Given two different auxiliary states representing the same kinematical state (i.e. 
having the same number of nodes), we can find two ordered sequences of operators 
that will produce the same abstract sum of triangulations when acting on the 
respective states. This is because the coordinates of the nodes just serve to label 
the links and their exact value is irrelevant for the geometrical interpretation of the 
action of Ui, a (x). But as noted above, in the limit that I — > 0, any triangulation 
that can be constructed by an ordered sequence of evolutions acting on a given 
auxiliary state, will be present in the expansion of P/. a . Hence, in this limit, the 
amplitude becomes independent of the choice of representative for the state \m). 
For the same reason the amplitude also becomes independent of the metric chosen 
in the regularisation of Pi >a 

To conclude we note that the triangulations constructed with the above defini- 
tions are unique, i.e. two different sequences of operators Ui ia [x\) ■ ■ ■ Ui, a {Xs) acting 
on the same auxiliary state will give rise to different triangulations. This follows 
essentially from the ordering of the operators as shown in the appendix. However it 
is also clear that not all triangulations are produced c.f. figure [|. This is a feature 
of the choice of definition of Ui, a {x). Indeed, it is possible to choose a different 
action of Ui. a (x) in which a(xi) corresponds to a gluing of a triangle to the links 
Xi and Xi+i. In this case all triangulations are produced albeit not uniquely. It 
remains to be seen if an action can be chosen such that we have both properties. In 
this respect it will be important to check to which universality class the continuum 
limit of the triangulations model constructed above corresponds to. As we discuss 
below the derivation of the Hamiltonian corresponding to this model suggests that 
we are summing over geometries that suppress the creation of baby universes. 



RELATING COVARIANT AND CANONICAL APPROACHES. 



13 




Figure 4. This triangulation cannot be constructed with an or- 
dered sequence of Ui >€ {x). This is because the top triangle has to 
be constructed after the bottom one. This cannot happen since it 
is situated to the right. In other words we would need a sequence 
Ui t a{xi)Ui^ a {xs), which is not allowed since £4 > X3. 



3.3. The Hamiltonian. One of the benefits of working in 1+1 dimensions is that 
we can use the definition of Ui ia (x) to derive the form of the continuum, aver- 
aged, Hamiltonian constraint, thus completing the link between the covariant and 
canonical pictures^. 

The Hamiltonian and U^ a (x) are related via: U l>a (x) = e - alHi ~ x \ In the following 
we will take the a and / limits simultaneously and hence set a — I. Thus for small 
a we have: U a (x) ps 1 — a 2 H(x). Let us define the averaged Hamiltonian as: 

H a = { dxH a (x) (25) 
Js 1 

Then 

aS a := / dxU a {x) » L^S 1 ) - a 2 H a , (26) 
Js 1 

where L(S 1 ) denotes the length of the spatial slice S 1 . The point is that the 
integrated evolution operator S a has a simple expression. It is the sum of the 
elementary evolution actions (adding of a triangle in two possible ways) on every 
node of a state \n). More specifically, the matrix elements of S a are given by 

(n\S a \m) = ye~ a2A ((5 m ,„_ 1 + S m<n+1 ) (27) 

The length ^(S 1 ) also has a natural operator action. Since the states \n) describe 
geometries of length na it follows that \n) are volume eigenstates, i.e.: 

V{Y)\n) = na\n). (28) 

Hence from equations (p6|), ( p7| ) and (|2^) we learn that the matrix elements of the 
discretised averaged Hamiltonian operator should be given by: 

(n\H a \m) fa — (5 m n - -e- a2A (5 m! „_ x + 5 m . n+1 )) (29) 
a 2 

The states \n) describe geometries of length N = an. We want to take the limit 
a — > while keeping TV fixed. In this way we obtain a continuum Hilbert space 
spanned by the states \N). We are looking for the continuum expression H of H a 



We thank Jan Ambj0rn for this observation. 
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which acts on this space i.e. we are interested in: 

Hf(N) = / dM(N\H\M)f(M) (30) 



Since the discrete states |n) = \N/a), |m) = \M/a) become increasingly better 
approximations to the continuum states \N) and \M) as a — > 0, we can replace the 
above by: 

Hf(N) = lim J2(n\H a \m)f(ma) (31) 



Using the fact that n = TV/a and m = M/a together with equation ( |29| ) we obtain 
l?/(iV) = Urn J(/(iV) - \e- a2 \f{N - a) + f(N + a))) (32) 

Expanding /(TV — a), f(N + a) to second order in a and also expanding the expo- 
nential this gives 

Hf(N) = hm^(/(iV)-i(l- a 2 A)(2/(7V)+a 2 -^/(7V)+0( a 3 )))(33) 

= N\f(N)-N±^f(N) (34) 
Hence we deduce 

H = XN-N l -±- 2 . (35) 

This is the same expression (up to a factor of 2) that has been obtained in the study 
of dynamical triangulation models [Q . It corresponds to a choice of triangulation 
models in which baby universe creation is not allowed^]. 

3.4. Discussion. We have succeeded in deriving the form of the averaged Hamil- 
tonian constraint operator from the regularised path integral for 1+1 quantum 
gravity using our projector expansion. This hinged crucially on the identification 
of the elementary evolutions with the addition of triangles. Moreover, this allowed 
us to naturally identify the small proper time step a occurring in the regularisa- 
tion of the projector with the edge length of the triangles. This gives rise to an 
interpretation to the integral over proper times [dX^x)] as a sum over triangulations. 

These features are not present in the original regularisation of the projector 
proposed by Reisenberger and Rovelli ^0|. They suggested an expansion of the 
exponential in the projector 

P=[[dT]jr (36 ) 

which as we will discuss later this was motivated by the form of the Hamiltonian 
constraint in loop quantum gravity. In order for this expression to be interpretable 
as a sum over triangulations it turns out that we need to identify H[T] with the 



7 Durhuus and Lee have shown || that lim a ^ (N/a\e- THa \M/a) = G(N, M;T). Here 
G(M, N;T) denotes a sum over sliced triangulations with T slices as depicted in figure [io| 
(see Jl| for details). From our perspective it is clear why this should be so. We have 
e~ TH = exp[— T J s x dxHa(x)] = [Il^gs 1 U a (x)] T . Hence the action of this operator corresponds 
to the stacking of T layers of triangles to a boundary triangulation. 
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integrated move operator S a . In this way multiple action of the Hamiltonian will 
produce all possible triangulations. 

This leaves the interpretation of the [dT(x)] integral unclear since a sum over 
all triangulations is already generated for every function T(x). Intuitively, the 
proper time should be reflected in the triangulation. Indeed in the dynamical 
triangulations approach it is clear that the proper time corresponds to the number 
of edges of triangles along a particular path. 

Additionally, a problem arises because there are many ways in which a given 
triangulation can be constructed by a sequence of stacking of triangles. All these 
possibilities will be included in the sum generated by the above expansion^]. This 
number is a non-trivial function of the triangulation and does not depend solely on 
t, the number of triangles in the triangulation. Hence the sum over triangulations 
will contain non-trivial symmetry factors that do not occur in the definition of 
G a (A,A'). 

Crucially, however, because one has made the identification H[T] = S a it will no 
longer be possible to derive the correct expression for the continuum Hamiltonian. 

3.5. Generalisations. It is straightforward to generalise our discussion to higher 
dimensions. Space-time triangles will be replaced by (n+l)-simplices, where (n+1) 
is the space-time dimension. These can be glued in (/ — 1) = (n + 1) ways to a 
given boundary triangulation, where / is the number of faces of the simplex. These 
moves will be identified with the action of the local evolution operator and in this 
way we can again interpret the projector as a sum over triangulations. This is 
described in more detail for (2+l)d in the next section. 

However in higher dimensions boundary states no longer only depend on the 
number of boundary simplices but also their connectivity, which is related in a 
complicated way to the limiting continuum geometry of a discrete state. This 
makes it very difficult to derive the action of a Hamiltonian differential operator. 
These issues are currently under investigation 

4. Spin-foam models and loop quantum gravity 

We will now discuss spin-foam models of BF theory and quantum gravity. In 
general, a spin-foam is a 2-dimensional complex with faces labelled by represen- 
tations and edges labelled by intertwining operators. These structures were first 
used in path-integrals for quantum gravity by Reisenberger | jl9| and now provide a 
unifying framework for many approaches to quantum gravity (see and |j| for 
reviews and references). 

In order to make contact with our discussion of dynamical triangulations, we will 
make use of the fact that the dual of a spin-foam is a triangulation of space-time 
with labelled edges. We will see how one can provide a definition of the projector 
(not sum over geometries) in terms of a sum over such space-time triangulations. 
In contrast to the dynamical triangulations framework, however, we also sum over 
the labellings of the edges in this approach. We will show how a new regularisation 
of the projector will lead naturally to this class of models. The focus will be on the 
(2+l)d case, and we indicate at the end how the results can be generalised. 

In 2+1 dimensions general relativity is equivalent to a topological BF theory 
due to the absence of local degrees of freedom or gravitons. In the Euclidean case 



This is avoided in our approach because the triangles are stacked in ordered sequences. 
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the fundamental variables are an su(2) Lie algebra-valued space-time connection A 
and an su(2)-valued one form B. The action is given by: 



where F(A.) is the curvature of A. By restricting to space-time manifolds of the 
form Ai = S x 1, where E is a compact 2-dimensional spatial manifold, we can 
pass to a Hamiltonian formulation. The configuration space variable is given by 
the restriction A of the connection A to E. The Hamiltonian is a sum of two 
constraints^. The first imposes St/(2)-gauge invariance and the second, also referred 
to as the curvature constraint, implies that F(A) — 0, i.e. that the connection is 
flat. A good overview and a comprehensive list of references can be found in [^|. 

It can be shown that the above is equivalent to a Chern-Simons gauge theory 
which was first used by Witten in |S3 to obtain a quantisation of 2+1 gravity. Here 
we will be using the discrete Ponzano-Regge |1| or Tuarcv-Viro Q approaches to 
quantisation as they will resemble the dynamical triangulations models discussed 
above. 

4.1. Quantum theory. The kinematical framework for the Ponzano-Regge model 
can be described by loop quantum gravity as first noticed by Rovelli ||22| . The 
kinematical Hilbert space Tikin is constructed from a space of gauge-invariant wave- 
functions on the configuration space of connections. This space is spanned by the 
so-called spin- networks [M, which are the natural generalisations of the Wilson 
loop. 

Spin-networks are functions of connections with support on graphs embedded in 
space, which in 2+1 we choose to be trivalent. More precisely, given any closed, 
oriented, trivalent graph P with n edges embedded in E and a labelling of the 
edges by representations of SU(2) denoted by j = {ji, ■-,jn}, we can construct 
a spin-network function^] ^r,j{A) as follows. Let us denote the holonomy of the 
connection A along the edge e S L by Z7 e [A]. This can be thought of as an SU(2) 
group element. Thus for every edge e labelled by a spin j e we can construct a matrix 
Pj c {U e [A])1 jl where pj e is the spin-j e representation of SU(2). We can associate the 
index a with the vertex to which the edge e is directed and the index b to the vertex 
at the other end of e. Taking the tensor product of these matrices gives us a tensor 
with three indices for every vertex in T. We can contract these indices with the 
unique invariant tensor in the tensor product of the representations on the edges 
incident on the vertex. This gives us the gauge- invariant function ^r,j{A). 

An inner product can now be imposed on the space generated by these functions 
by demanding that spin- networks be orthonormal, i.e.: 



Completion with respect to this inner product gives us the kinematical Hilbert 
space Hkin- 



9 Note that there is not need for a separate 2-diffeomorphism constraint as such transformations 
are already generated by the curvature constraint. 

10 Note that because we have chosen graphs to be trivalent we do not require a choice of 
intertwiner for the vertices. In 2+1 this does not imply a loss of generality since the theory can be 
solved exactly in terms of these spin-networks. Whether this is also possible in 3+1 dimensions 
with 4 valent graphs is unclear. 




(37) 



(38) 
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We can also describe the trivalent spin-networks in terms of their dual trian- 
gulations A, making the link to earlier discussions. Each edge of A is labelled by 
the spin on the edge in T intersecting it. We will switch between these two per- 
spectives whenever convenient and spin-networks will be labelled either by graphs, 
triangulations or not at all. 

To proceed we now need to implement the Hamiltonian or curvature constraint 
H and reduce to the final, physical, Hilbert space H p h ys - As before this is done by 
constructing the generalised projector. In contrast to the previous section one can 
give an exact expression of this projector in terms of sums over triangulations. 

As shown by Ooguri [jl6|, the Ponzano-Regge model defines the projector for 
2+1 quantum gravity. This is most easily introduced in terms of triangulations ▲ 
of M. = £ X [0, 1]. We will distinguish between the interior edges e and the edges 
e of the boundaries Ao and Ai of ▲. All edges are labelled by representations j e 
of SU(2). Hence we can associate to every tetrahedron T c A a normalised 6- J 
symbol denoted by T{jr}, where jt is the set of spins on the boundary edges of 
T, which can include edges that lie on the boundary of A. To every triangulation 
▲ with given boundary Ao, Ai wc can now define an amplitude Z A (Ao,Ai) as a 
sum over all possible labellings of the interior edges that are compatible with the 
labelling of the boundary: 

z A (A 0> Ai) :=nv / 2?m^n( 2 ^+ i )n r {^}- ( 39 ) 

e je e T 

This sum can be infinite and needs to be regulated. This can be done by intro- 



ducing a cut-off as is done by Ooguri in 1 16 by using quantum groups. The latter 
approach gives the Tuarev-Viro model |33 | which is related to quantum gravity with 
a cosmological constant. This will not concern us in the following and we assume 
that some choice of regularisation has been made to make the above expression 
well-defined. 

Given two spin-networks \&n and ^Pi, dual to triangulations Ao and Ai, we can 
now define the physical inner product as a sum over triangulations A that have 
boundary Ao + Ai: 

(*l,P*0>«n = £)^(Ao,Ai) (40) 

A 

Using identities between the 6-J symbols one can show that Za(Ao, Ai) is in fact 
independent of the triangulation ▲ and depends only on Ao and Ai. Hence the 
matrix elements of P are given by oo x Z(Aq, Ai), where Z is evaluated using any 
chosen triangulation with the correct boundary data. Such an infinite factor is a 
common feature in the RAQ procedure whenever the constraints do not form the 
Lie-algebra of a compact group^j |J . The factor can be dropped by redefining the 
physical inner product and we have the final resultp]: 

^o)phys = Z(Aq, Ai). (41) 



11 Recall that the projector averages over the action generated by the Hamiltonian constraint. 

12 The reason that we chose to include a sum over triangulations initially is that it will play 
an important role in the following. The above procedure is also necessary for gravity in higher 
dimensions where we no longer have triangulation independence. 
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A crucial property for the identification of the amplitude Z with a projector is 
the property: 

Z{hx, A J )Z(A J , A ) = Z(Ai, Aq), (42) 

j 

where the sum is over all colourings Aj of an intermediate triangulation A. 

As has been shown by Ooguri |L6]], this definition of P leads to a well-defined 
theory of quantum gravity, which is equivalent to Witten's Chern-Simons approach. 
If we formulate the above in terms of spin-networks instead of their dual triangu- 
lations we obtain the spin- foam model of 2+1 loop quantum gravity. A spin- foam 
is simply the 2-skeleton of the dual 2-complex of a given triangulation. 

The sum over spins on a triangulation can be given an interpretation as a sum 
over geometries or path integral in several ways. Originally the Ponzano-Regge 
model was an ad hoc modification of the Regge model of gravity. Here one de- 
scribes (piecewise linear) metrics by an assignment of edge lengths and angles to 
the tetrahedra in a triangulation of space-time. By restricting edge lengths to take 
on only integer and half-integer values we obtain the Ponzano-Regge model. 

The actions of the two models can be related in a semi-classical limit where we 
take spins j to be large. In this limit the amplitude for a single tetrahedron given 
by the 6-J symbol can expressed as: 

T{j T } ~ (6^Vr)" 1/2 (e- 4(Si?+ " /4) + e i{SR+ *!^), (43) 

where Vt is the volume of the tetrahedron sides of length L — ^(j + |) and Sr 
is the Regge action. Hence it is natural to identify spins with geodesic lengths 
in the Ponzano-Regge model. The fact that the 6-J symbol describes a sum of 
two exponentials with opposite signs is related to the fact that the Ponzano-Regge 
model defines a true projector and hence not a sum over geometries as we comment 
on again below. 

A perhaps better understanding of the origin of the discreteness in the edge 
lengths comes via the spin foam approach. It is possible to show that Ponzano- 
Regge model arises if one discretises the path integral: 

j ' [dB][dA]e- tS ^' A \ (44) 

where the action is defined in equation (^) . The spins labelling edges then arise af- 
ter an integration over B, as the Fourier components of the connection A restricted 
to the edges of the triangulation. Furthermore, one can construct a length operator 
for loop quantum gravity and one finds that the spins are precisely proportional to 
the quantised edge lengths as first noticed by Rovelli in p2f . 

Since spin-foam models involve a sum over edge lengths, a continuum limit can- 
not be taken in the same way as in the dynamical triangulations models where 
edge lengths are taken to 0. The implications of this will be discussed below. It 
turns out that we need to take the scale of the edge lengths to 0. Note however 
that because the special cases of BF theory and 2+1 quantum gravity are topolog- 
ical these models can be defined exactly in terms of a single triangulation without 
taking a continuum limit as described above. 

Note also that the sum over geometries that is described by the Ponzano-Regge 
and more general spin-foam models is not gauge fixed. In other words different spin 
assignments of a triangulation can correspond to the same geometry. The rationale 
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Figure 5. A boundary triangulation Ao- 




Figure 6. Gluing one face a tetrahedron to Ao- Only the relevant 
spins are labelled. 

behind the spin foam approach is that instead of a picking a representative of each 
gauge equivalence class and integrating only over these wc arc in effect averaging 
over the equivalence classes. This will be important later. 

4.2. The projector and elementary moves. In order to derive the sum over 
triangulations picture from the formal definition of the projector we now describe 
how one can construct Z^(Aq, Ai) for any ▲ and Ai given an initial triangulation 
Ao by an iterative sequence of elementary moves. These are the Pachner moves in 
2+1 dimensions by which all triangulations can be related. 

There are 3 ways of gluing a tetrahedron T onto a given triangulation Ao of S. 
These are shown in figures |[ and H: 

1. We can attach one face of T to a triangle in Ao- To this tetrahedron we 
associate the factor 



where we have written out the 6- J symbol and k,l,m are the spins on the 
edges of T that are not glued to Ao 
2. We can attach two faces of T to Ao- This is associated with a factor 



where k labels the edge that lies in the interior of the two triangles glued to 
Ao and I is the spin on the opposite edge. 
3. We can glue three faces of T to Ao- In this case we have a factor of 



where k,l,m label the edges in the interior of the three triangles attached to 




(45) 




(46) 




(47) 



A . 
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Figure 7. Gluing two faces of a tetrahedron to Aq. 




Figure 8. Gluing 3 faces of a tetrahedron to Aq- 



Once we have completed the gluingp^ we obtain a new boundary triangulation 
A'. Using this boundary we can repeat the above procedure. In this way we can 
construct any triangulation A with any boundary Ai given an original boundary 
Ao- By taking the product of the factors associated to the tetrahedra we recover 
the amplitude Z A (A , Ai). 

Equivalcntly such a procedure can be viewed as constructing a sequence of spin- 
networks via a sequence of operator actions, which can be formalised as follows. We 
introduce the 3 move operators mi, m 2 and 7713, that act locally on spin-networks: 



h 



mi 



Sfe A.m. T 



3 1 
I 



32 

m 



33 
k 



v /(2/c + l)(2; + l)(2m + l) 
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Jl js \ X J(2k + l)(2l + l) ■V J -<^ 
J4 32 k J v A J\ 



m 3 



'Jl J3 



T 



31 32 33 

I m k 



v /(2/c + l)(2; + l)(2m+l) . 



Jl 



By a finite sequence of such moves one can relate any two trivalent spin-networks. 
Note that mi, m 2 and 7713 are precisely dual to the three possible gluings of tetra- 



hedra as shown in figure 
equations (j45|), (Eq) and <k 



The amplitudes also correspond to those given in 



13 In terms of spin-foams we are gluing together elementary atoms as is described in ^jj. The 
different ways of gluing triangles correspond to the choice of which boundary edges of the atom 
we wish to glue to the given spin-network. 
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Figure 9. The action of the elementary moves on spin networks 
are dual to the addition of tetrahedra. Tetrahedra are depicted 
by the thin lines and are viewed from the top. For example, we 
can view top left figure as boundary triangulation consisting of 3 
triangles with its dual graph. By gluing a tetrahedron to all three 
triangles we obtain a new boundary triangulation consisting of just 
one triangle as shown on the right. This corresponds to the move 
fri3 acting on the dual graphs. 



Hence any coloured triangulation ▲ can be viewed as a given sequence of such 
elementary moves. 

As with the dynamical triangulations models, we wish to use the fact that the 
projector can be expressed as an ordered sequence of operations in order to derive 
the above action from the formal expression of the projector. This will clarify how 
the elementary move operators mi, rhi 1 rh^ are related to evolution operators and 
in consequence the Hamiltonian. 

4.3. Regularisation of the projector II. We wish to provide a regularisation of 
the generalised projector as a product of elementary local operators in analogy to 



the approach in section f^lj. These elementary operators are to be identified with the 
moves rhi , to 2 and to 3 . But note that in contrast to the dynamical triangulations 
case the these moves do not involve a fixed proper time, instead the describe a sum 
over all possible proper times. This needs to be reflected in our expansion. 

It turns out that the correct starting point for our calculation is the following 
expression for the projector: 

P:= j[dN{x,t)]e- l ^ dt S d 2 xN(x,t)H(x)^ ( 4 g) 



where we are integrating over all lapse functions on the cylinder £ x [0, 1]. The 
reason we need to extend the range of integration as compared to our previous 
definition in equation (||) is that the path integral that we wish to recover is not 
gauge fixed as mentioned already above. This is reflected in the fact that the two 
definitions of the projector are related by: 

[dN(x,t)]e- i ^ dt ^ xN(x ' t) ^ x) =Af [[dT(x)]e- iA W, (49) 
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where T(x) = /„ dtN(x, t) and TV is an infinite normalisation factor, which should 
correspond to the volume of the gauge equivalence classes. Finally, we note that in 
this section we integrate over the entire range of lapse functions since we want to 
describe a genuine projector. 

Let us now proceed with the regularisation of P by replacing integrals by Rie- 
mann sums. First we consider the integral over the coordinate time t and we replace 
the integration over functions N(x, t) by a product of integrations over spatial func- 
tions N(x). This leads to the following regularised expression: 

[dN(x, i)]e- 4 £to £ 1 d ^ N i^H{x) ( 50 ) 



= JJ f{dN(x)}e- i ^ d2xN ^" (x \ (51) 
t=o J 

where k = 1/e. Now we proceed analogously with the remaining spatial integral 
in the exponential. To do this we need to introduce a square lattice with a set of 
vertices V as before. 

P ltt = Yl f[dN(x)]e- i€ ^-evi 2 N(x)H(x) (52) 

+ — n •* 

dTe - t el*TH( X ) (53) 



t=0" 
fe 



nn 

t=0x£V 
fe 



— oo 



roc 

= II II / dT (e- id2TH(x) + e id " rB ^) , (54) 
where we are left with a regular integral over TeR. Let us now set 

/>oc 

U lte (x) := / dT (e~ id2TH W + e ^ TH ^) (55) 



then we have the final expression: 

k 

he = n n uuz) ( 56 ) 

t=0x£V 

This expression is similar to equation (^) and can be interpreted as a sum 
over triangulations. There are however essential differences, most notably in the 
definition of Ui tt . Let us take these in turn. 

4.4. Interpretation. The operator U^ e (x) we have defined above is given by a 
sum over local evolution operators for all proper times. Hence it is natural to make 
the identification: 

£/j, e (x) — mi+m 2 + m 3 (57) 

The argument x of the operator Ui^(x) determines the location of the action of the 
move operators rhi on a given spin network. The integral over T in the definition 
of Ui tt essentially reflects the sum over spins in the definition of the elementary 
moves. This identification is further supported by the fact that Ui :<L contains two 
terms involving an evolution in opposite proper time directions. This mirrors the 
semi-classical interpretation of the 6-J symbol given a sum of two exponentials of 
the Regge action described in equation (Ba). Here we can see that this is a result 
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Figure 10. A (l+l)d triangulation with a slicing structure. 

of the integration range in the integral over lapse functions in the definition of the 
projector and hence a consequence of the fact that we are not dealing with a sum 
over geometries but a genuine projection operator. 

Once we have made the identification of Ui t6 with the elementary moves we 
can see how equation (|5^) describes a sum over triangulations. Each operator 
Ui t e(x) adds a tetrahedron in all possible ways and with all possible spins to a 
boundary triangulation at location x. As we take the limit I — > the product over 
positions x &V then describes the addition of tetrahedra to every position on the 
initial boundary, such that we obtain a single layer of tetrahedra. We then repeat 
this at the next co-ordinate time step. The result is that we obtain a sum over 
triangulations which each have a layered structure as depicted in figure [To| . 

These layers can be interpreted as the slices {t} x S in the foliation £ x [0, 1]. 
The spins on the edges of the tetrahedra describe the multi-fingered proper-time 
separation between the slices. The sum over spins then includes all possible proper- 
time separations between initial and final boundaries. 

We have so far not addressed how the parameter e is related to the moves rhi 
and how this continuum limit e — > is to be interpreted. Indeed, as mentioned 
above we can no longer identify e simply with a proper time separation since we 
are summing over all possible proper times for every tetrahedron. Instead we see 
from equation ( ]55| ) that e multiplies the proper time T and hence is responsible for 
determining the scale of the sum over proper times. This indicates that e should be 
identified with the scale associated in the relation between spins and edge lengths of 
tetrahedra. Indeed, in a continuum limit we would like the edge lengths to become 
small in some units. Hence, we should set physical lengths proportional to spins 
via: L ~ ej. 

Remarkably this parameter arises naturally in 3+1 loop quantum gravity. It is 
the Immirzi parameter [ fi"l[ fj5| 7 , which describes a one-parameter family of clas- 
sical canonical transformations which in turn give rise to a one-parameter family of 
inequivalent quantum theories. This parameter relates spin labels of spin networks 
to the spectra of physical geometrical operators. It determines that the eigenvalues 
are given in units of {^lp) a where l p is the Planck length and the power a is the 
dimension of the geometric observable we are considering. Hence, we deduce that 
we should make the identification 7 = e and a continuum limit of loop quantum 
gravity is obtained by taking the Immirzi parameter to 0. This is especially inter- 
esting as it has been proposed in an entirely different context by Bojowald who 
studied reduced cosmological models using loop quantum gravity. 

We note again that BF theory and (2+l)d gravity are special in that they are 
topological and the limit e — > does not in fact need to be taken to define the 
theory. This can be interpreted as stating that it does not matter at which scale we 
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approximate a continuum geometry. Since there are only a finite number of degrees 
of freedom any discrete approximation will capture the full physical content of the 
theory. 

4.5. The Hamiltonian. In the previous section we have provided a well motivated 
expansion of the projector describing the Ponzano-Regge sum over triangulations. 
Some of the attractive features of this approach were: 

• The natural identification of the sum over spins as a sum over proper-times; 

• The explanation of the semi-classical limit of 6-J symbols in equation (^) as 
a result of the integration range in the generalised projector; 

• The identification of the continuum limit as a limit in which the scales of edge 
lengths are taken to 0; 

• The reflection of the fact that the Ponzano-Regge model does not describe 
a gauge fixed sum over metrics in the definition of the projector given in 
equation (f48|). 

In this section we discuss a further motivation for our approach, which is closely 
linked to questions concerning the Hamiltonian. The point is that the local move 
operators fhi leave physical states invariant, i.e.: 

(*, rhxfrjphys = rh2$) p hys = m 3 $) phys = (*, $) phya (58) 

for all "P, $ This is shown by Ooguri in |u|. Let us take to 2 as an example. We 
have: 

m 2 $A =^T{J T } v /(2fc + l)(2j + l)<I> Afc =^Z(A fe ,Ao)$ At , (59) 

k k 

Here $ Ao is the spin-network dual to the coloured triangulation Ao and A& is the 
triangulation which is obtained from Ao by adding a tetrahedron with external spin 
k as shown in figure |^. The last equality in the above equation follows from the 
definition of Z and the fact that, due to the triangulation invariance of the model, 
the transition amplitude between the states $ Ao an d 3>aj is given by the amplitude 
of a single tetrahedron, which in turn is the same as the amplitude of the to,2 move. 
Hence, it follows: 

(* Al ,m 2 $AoWs = J2 Z ( Al > Afc)Z(A fe , A ) (60) 

k 

= Z(Ai,A ) = <*,$W, (61) 

where we have used the projection property ( |42| ) of the amplitude Z. 

This can also be understood as follows. The constraints of 2+1 general relativity 
impose the condition that the connection is flat. Physical states can be though of 
spin-network functions evaluated on flat connections. But it can be shown that 
the action of the elementary moves on spin-network functions restricted to flat 
connections is the identity. 

Physical states are orbits of kinematical states under the exponentiated action 
of the Hamiltonian constraint. Hence the elementary moves map between states in 
Hkin that are related by such action and they can thus be viewed as elementary 
time-evolutions. Thus it is natural to identify the moves with the evolution operator 

Ul,e- 

Note, however, that this is in contrast to earlier assumptions on the action of the 
Hamiltonian in loop quantum gravity. Let us recall that in the Reisenberger-Rovelli 
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approach the exponential in the projector is expanded as follows: 

f [dT{x)]e-^ = [ [dT{x)\ f] t^H[T]\ (62) 
•* t=o 

The idea was to identify H [T] with the total move operator 

H[T] = S := J dx(a 1 m 1 + a 2 m 2 + a 3 m 3 ), (63) 

where a\,a 2 and a 3 are arbitrary coefficients that need to be fixed. This was 
motivated because the action of Thiemann's Hamiltonian ]32] | in 3+1 dimensions 
resembled the move rh\ albeit with a different amplitude. The other moves need 
to be included in the Hamiltonian if the multiple actions of the Hamiltonian are 
to give rise to all triangulations. This is because the integral over the scalar lapse 
functions does not alter the topological action of the operators. The search for such 
a Hamiltonian, also referred to as crossing symmetric, has met only with partial 
success [0. In particular it is not clear how the move m 2 , which is thought to 
be crucial for long range propagation of excitations can be derived from the 
quantisation of a classical expression of a Hamiltonian. 

But from the discussion above one can understand why this might be the wrong 
approach. Since the local moves leave physical states invariant it is impossible 
that any combination of them will annihilate physical states as would be required if 
they were to correspond to a Hamiltonian constraint. This is because the integrated 
move operator will be a sum of elementary operators acting at all possible locations 
in the triangulation underlying a spin-network. Each single action leaves physical 
states invariant. This implies the following for the action of the Hamiltonian defined 
in equation ( |63| ) on spin networks ^ and $ 

F[T]*[$] := / da;($|P(aimi + a 2 m 2 + a 3 ro 3 )|#) (64) 
= (ain x (*) + a 2 n 2 {^) + a 3 n 3 (*))^,$, (65) 



where we have used equation (58) and the orthogonality of spin- network states in 
the kinematical inner product. The coefficients m, n 2 , n 3 count the number of 
times that the operators mi, rh 2 and m 3 respectively can act on the state ^ . But 
since these coefficients are state dependent and the coefficients a, b, c are not, the 
action of H[T] cannot be made to vanish on all physical states. 

So from our perspective it seems clear that the difficulties in generalising the 
definition of the loop quantum gravity Hamiltonian stems from identifying the move 
operators with Hamiltonian actions and not evolutions. This is in line with views 
expressed by Markopoulou fl4| , where evolution of spin networks is described in 
terms of the Pachner moves listed above. In so far as these models should reproduce 
quantum gravity in 2+1 dimensions our work gives them a precise motivation in 
terms of canonical loop quantum gravity and shows that the amplitudes for these 
transitions should be given by the 6-J symbols as described above. 

The problem we face is to find a Hamiltonian which can reproduce the elemen- 
tary moves when exponentiated. This seems a difficult task. One possibility is that 
it would be possible to quantise the action generated by the classical Hamilton- 
ian directly, in the same way that it is only possible to solve the diffeomorphism 
constraint by quantising the action generated by it. 
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Another option that is currently being explored is to make use of the close con- 
nection to the classical formulation of 2+1 gravity in terms of 't Hoofts lattice 
model p^]. Here states are described by trivalent graphs with labellings corre- 
sponding to edge lengths. It can be shown that Hamiltonian evolution increases 
and decreases the edge lengths. Crucially, edge lengths can evolve to zero at which 
point a change in the topology of the graph can occur. These possible changes have 
been tabulated in p^j and, in the absence of matter, they turn out to be precisely 
the three elementary evolution moves. Hence it is possible that Hamiltonian action 
should act on edges by increasing and decreasing the values of the spins by adding 
loops of spin 1 to all closed loops in the graph of a given spin-network. 

4.6. Generalisations and outlook. Our results hold for BF theories in any di- 
mension, which are all topological field theories. These models arc also expressed in 
terms of a single triangulation, where the 6-J symbols are replaced by their higher 
dimensional analogues. 

The situation looks more complicated for quantum gravity since in space-time 
dimension higher than 3 the theory has an infinite number of degrees of freedom, 
and is no longer equivalent to BF theory and no longer topological. However there 
still exist spin-foam models that attempt to define a path integral for the theory 
both for Lorentzian and Euclidean signatures (see [[l7] for a review). In this case 
the transition amplitudes are no longer triangulation independent and we need to 
include the sum over triangulations to restore covariance. However, our results 
have shown that the sum over triangulations should be restricted to triangulations 
with a slicing structure. Such restricted sums have been studied in the context of 
Lorentzian dynamical triangulation models Q and there is mounting evidence that 
they have better behaviour than sums over generic triangulations. 

In addition, as soon as spin-foam models are no longer topological, taking a 
continuum limit will become important. Our work indicates that in loop quantum 
gravity and related spin-foam models this limit involves considering large spins 
while taking the Immirzi parameter 7 to 0. This has been suggested independently 
by Bojowald in the context of cosmological models [Q. In addition, the study of 
the geometric operators area, volume and angle around a vertex of a spin-network 
carried out by Major and Seifert Jl3| has revealed that in order to approximate 
continuum geometries resembling what we observe, we need to consider a large 
spin limit. 

We note also that for non-topological theories we can no longer argue that the 
elementary moves leave physical states invariant. Presumably again this is linked 
to the need to take the continuum limit in defining the evolution operators. The 
physical consequences of this deserves to be investigated in more detail. 
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Appendix 

We show that if a triangulation is generated by the action of an ordered product 
Ui } a{xi) ■ ■ ■ Ui^ a {xs) acting on an auxiliary state, then this is the unique sequence 
of operators that will produce this triangulation, where Ui yCL {x) is defined in equa- 
tion ( p2[ ) . We proceed by induction over the number of triangles in a triangulation. 

Clearly the addition of one triangle to a boundary triangulation corresponds to 
a unique action of on Ui ja (x). The position x is given by the co-ordinate of the 
node on the right edge of the triangle glued to the boundary. 

Let us assume that any triangulation with n triangles has a unique sequence of 
evolution operators corresponding to it. Now let us consider a triangulation with 
n + 1 triangles, which can be described by at least one sequence of operators. This 
triangulation will have a set of triangles that is attached to the boundary. This set 
inherits an ordering from the order of the coordinates of the nodes on the edges 
that the triangles are glued to. Hence, the triangulation will contain a unique 
triangle, which is defined by the fact that it is attached to a boundary edge with 
the largest value of the node coordinate. Since the products of evolution operators 
in the expansion of P^ a are ordered this triangle has to be constructed first. Hence 
the first operator in any sequence of operators corresponding to our triangulation 
is fixed and unique. 

But now we need to add only n more triangles to obtain our triangulation. By 
assumption there is only one sequence of operators that will describe this triangu- 
lation with n triangles. 

Hence we conclude that every triangulation with n + 1 triangles has a unique 
corresponding product of operators Ui, a {xi) ■ ■ ■ Ul,a( x n+i)- qed. 
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